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ABSTRACT 

We hereby report on a sensitive search for radio continuum observations from a sample of 34 
Mira and semi-regular variable stars. The main aim of this survey was to search for thermal 
free-free emission from post-shock ionised gas. Thirty-four stars were observed at 3- and6-cm 
using the Australia Telescope Compact Array. Radio continuum emission was detected from 
one source only, the symbiotic Mira R Aqr. No continuum emission was detected from the 
other sources, with three-sigma upper limits of typically 0.3 mJy. From the upper limits to the 
radio flux densities, we have found upper limits to the gas brightness temperatures near two 
stellar radii at a characteristic size of 5 x 10 13 cm. Upper limits to shock velocities have been 
estimated using a shock model. For the 1 1 nearest sources in our sample we obtain brightness 
temperatures below 6 000 K and shock velocities below 13 km s _1 . For 1 1 out of 14 sources 
with previously published detections, the radio brightness temperatures are below 4 000 K. 
For an upper limit of 4 000 K, we estimate that the shock velocities at two stellar radii are 
below 10 kms -1 . 

Key words: circumstellar matter - infrared: stars. 



> . 1 INTRODUCTION 



Mira variables are optically- visible long-period variable stars with 
large-amplitude pulsations (Ami, > 2.5'") and stellar periods over 
100 days. The stars have reached an advanced stage of evolution 
on the asymptotic giant branch (AGB). Miras typically have photo- 
spheric radii of ~ 350 R©, stellar luminosities of 10 3 -10 4 L© and 
effective temperatures of ~2 500 K. 

During their evolution on the AGB, Miras develop increas- 
ingly strong stellar winds, in which mass is lost from their outer 
hydrogen envelopes, with mass-loss rates M above 10 -7 M© yr _1 
and outflow velocities of typically 8kms _1 . Such mass-loss plays 
a critical role in their subsequent stellar evolution, leading eventu- 
ally to the formation of a white dwarf star and planetary nebula. 

The circumstellar mass-loss of Miras and other related AGB 
stars is driven by radiation pressure on solid grains, which form 
in the extended hydrogen atmospheres at heights of several stel- 
lar radii. Once formed, the grains are rapidly accelerated out- 
wards, while grain-gas collisions lead to a n expanding circumstel- 
lar envelope dGoldreich & ScoviHe|[l976l) . The radiation pressure 
model is well established but the physical conditions and mass- 
loss processes within the grain formation zone, are poorly under- 
stood. It is generally accepted that the high mass-loss rates are re- 
lated to the stellar pulsations and strong shock waves which prop- 
agate outwards during each pulsation cycle, increasing the mass 
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density in the upper s tellar atmospheres to a height wh e re dust 
formation occurs (e.g. IWood & KahiJ|l977l : IWoodl 1 19791 : iMorrij 
1 19871 : lAnandarao. Pott asch & V aidval 1 19931) . Shock velocities in 
the optically visible photospheres have been estimated from hy- 
drogen B aimer a nd metallic emission lines to be between 10 and 
90 kms" 1 (e.g iFox Wood & Dopital 1 1984 iGillet et al.l Il985bl : 
iRichter & WoodkOOll) . 

iRudnitskiil Jl99Ct> and iLuttermoser & Brown! \ 19921) argued 
that at centimetre wavelengths, a post-shock partially-ionised layer, 
extending to several stellar radii, would emit optically thick ther- 
mal free-free emission. For a Mira variable at a distance of sev- 
eral hundred parsec, the expected peak flux density from such a 
layer would be several milli-Janskys which is easily detectable. In 
addition, thermal emission from the stellar photosphere and non- 
thermal radio emission associated with stellar flare activity might 
also be detectable. 

So far, radio continuum emission has been detected from 
approximately 1 5 Mir as and semi-regular variables (Section 4). 
iReid & Mentenl i ll 9971) reported the detection of optically-thick 
thermal continuum emission from six nearby long-period variable 
stars. Their VLA radio flux densities are a factor of about two larger 
than e xpected for emission fr om the optical stellar photospheres 
alone. Rei d & Mentenl i 19970 provided a convincing model for 
their data in which the radio emission occurs from a 'radio photo- 
sphere', of about twice the optical diameter, where the gas tempera- 
ture and density are approximately 1 600 K and 1.5 x 10 12 cm -3 re- 
spectively. For these sources, the low radio brightness temperatures 
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are not consistent with strong shock activity in the radio emitting 
regions. However, earlier radio stud ies of several sources, includ- 
ing R Aql, o Cet and CW Leo (eg. I Woods worth & Hughesll 19731 
Spergel, Giuliani & Knapp 1983; Estalella, Paredes & Rius 
Knapp et all 1 19951) reported much stronger radio flux den- 



1977 




1983 





sities, consistent with higher brightness temperatures. 

In this paper we describe a survey in which we used the Aus- 
tralia Telescope Compact Array (ATCA) to search for radio contin- 
uum emission at A = 3 and 6 cm from a sample of 34 Miras and 
semi-regular variables. The primary aim of this study was to use 
radio continuum flux densities or their upper limits as a means of 
investigating whether high radio brightness temperatures, consis- 
tent with strong shock activity, may be present in some sources. 



2 SOURCE SELECTION 

2.1 Selection Criteria 

A sample of 34 Miras and semi-regular variables was selected us- 
ing the dat a in the Fourth edition of the General Catalogue of Vari- 
able Stars dKholopov etaTIl 1986-1989h . Of the 34 sources, 26 are 
southern Mira variables with declinations below 0° and visual mag- 
nitudes at maximum brightness of m v < 7 m . One of the selected 
Miras, R Aqr, is a well-known symbiotic binary star. In addition, 
three northern Miras, U Ori, U Her and R Aql were also observed: 
U Ori and U Her as they have well-studied OH maser pro perties 
dChapman & Cohenl 1 19851 : IChapman. Cohen & Saikia]|l99ll) and 
R Aql as it has been previously detected in radio continuum emis- 
sion. 

The sample also included five southern sources, R Pic, L2 Pup, 

V Hya, W Hya and S Pav, which are classified as semi-regular 
variables but have well-defined light variations. The semi-regular 
variables have smaller amplitude pulsations and lower mass-loss 
rates than the Miras and might be expected to have weaker shocks. 

V Hya and W Hya, however, have previou s ly been detected i n radi o 
continuum dLuttermoser & BrowrJ 1 19921 : iReid & Menter] 1 19971) . 
L2 Pup is one of the closest long-period variable stars, with a 
maximum- visual brightness of 2.6 m . 

Table 1 lists properties of the observed sources. The columns 
of the table are 

(1) source name 

(2) IRAS identification 

(3) variability type 

(4) P, the stellar period 

(5) (m v ), the mean visual maximum 

(6) D, the stellar distance 

(7) <p, the phase of the visual light curve 

(8) ft.so, the 3<r upper limits to the 4.8 GHz (6 cm) flux density 
(Section 3) 

(9) 5*8.64, the 3a upper limits to the 8.64 GHz (3 cm) flux den- 
sity (Section 3) 

(10) Tt, the upper limit to the brightness temperature of the ra- 
dio emission, for an adopted radio size of 5 x 10 13 cm (Section 
4.2) 

(11) v s , the upper limit to the shock velocity. 

2.2 Stellar Distances 

In Table 1, the stellar periods, mean visual maxima and phases 
have been taken from optical light curves provided by the Amer- 
ican Association of Variable Star Observers (AAVSO). For 10 of 



the nearest sources (o Cet, R Hor, R Lep, L2 Pup, R Car, W Hya, 
R Nor, R Aql, RT Sgr, and R Aqr), the stellar distances have 
been determined directly using trigonometric parallax measure- 
ments recently provided by the Hipparcos astrometry mission. For 
these sources, the distance errors given in Table 1 are determined 
from the standard errors in the stellar parallaxes as given in the 
Hippa rcos Catalog's] dl997l) . For one nearby source, R Hya, Hippar- 
cos parallax data were not available. 

For the more distant sources, with distances above ~320 pc, 
the Hipparcos parallax measurements are too small to provide re- 
liable distances. For 23 sources, the distances have been calcu- 
lated using the relationship between stellar per iod and absolute 
visual magnitude at mean maximum gi v en by IClavton & Feastl 
dl969h : IFqv, Heck & Mennessier M. -O. I dl975h . Interstellar ex- 
tinction corrections were c alculated iterat i vely using Van Herk's 
model of visual absorption dVan Herkll965l : IClavton & Feastl 19691) 
with: 



A v = 0.14| csc6|[l - exp(-0.01D| sin6|)], 



(1) 



where b is the galactic latitude, D is the source distance in pc and 
Ay is the visual absorption in magnitudes. The stellar distances 
were then calculated in the standard way with 



51ogD = (m v ) - M v 



(2) 



where my- and Mv are the apparent and absolute V-magnitudes 
at mean maximum. Mira distances calculated in this way are accu- 
rate to typically 30% with the uncertainties arising largely from the 
intrinsic scatter in the period-luminosity relation. Because of their 
irregular light variations, larger distance errors are expected for the 
semi-regular variables. For one source, V Hya, the distance has 
been determined using the near-infrared I-f -band period-luminosity 
relation for oxygen-rich Miras in the LMC, given by Feast et al. 
(1989), with a distance modulus t o the LMC of 18.55 and averag e 
value of Jf-magnitudes given by iNeugebauer & Leightonl dl969h : 
iGezari, Schmitz & Meaddl984h : lFouque et al.l dl992l) . 

For nine of the ten nearest sources, a comparison of the dis- 
tances obtained from the Hipparcos data with distances determined 
from equations (1) and (2) shows good agreement within the likely 
errors. For one source, RT Sgr, the Hipparcos distance of 133 pc is 
considerably smaller than the distance of 408 pc determined from 
the P—L relation. The reason for this discrepancy is not clear. 



3 OBSERVATIONS AND RESULTS 

The observations were taken on 1995 November 3-5, using the 
Australia Telescope Compact Array (ATCA) with a maximum base- 
line of 6 km. The ATCA is a synthesis instrument consisting of six 
22-m antennas on an east-west track, located near Narrabri in New 
South Wales. The selected sources were simultaneously observed 
in two frequency bands centered at 8.64 GHz (3 cm) and 4.80 GHz 
(6 cm) with a bandwidth in each case of 128 MHz. Each source was 
observed for a total time of approximately 75 minutes, split into 
five scans of 15 minutes with the scans distributed over the sidereal 
time range of the source. The data were corrected for atmospheric 
amplitude and phase variations using observations of nearby strong 
continuum sources, and the absolute flux density scale was cali- 
brated using the primary calibrator source 1934—638, which was 
taken to have 6 and 3 cm flux densities of 5.83 and 2.84 Jy, respec- 
tively. 

The data reduction was carried out in a standard way using the 
AIPS software package. After editing and calibrating the data, the 
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Table 1. Stellar Properties, Upper Limits to Radio Fluxes and Brightness Temperatures 

Upper limits 

Star IRAS Type P (m v ) D* ip S4.8O ^64 T** ~ s 

(days) (mag) (pc) (mJy/Beam) (K) (kms _1 )t 



S Scl 


00128- 


-3219 


M 


Cet 


02168- 


-0312 


M 


UCet 


02313- 


-1322 


M 


RHor 


02522- 


-5005 


M 


RRet 


04330- 


-6307 


M 


RCae 


04387- 


-3819 


M 


RPic 


04448- 


-4920 


SR 


RLep 


04573- 


-1452 


M 


UOri 


05528+2010 


M 


VMon 


06202- 


-0210 


M 


L 2 Pup 


07120- 


-4433 


SR 


RCar 


09309- 


-6234 


M 


SCar 


10077- 


-6118 


M 


VHya 


10491- 


-2059 


SR 


XCen 


11466- 


-4128 


M 


UCen 


12307- 


-5422 


M 


UOct 


13182- 


-8357 


M 


RHya 


13269- 


-2301 


M 


WHya 


13462- 


-2807 


SR 


RCen 


14129- 


-5940 


M 


RSLib 


15214- 


-2244 


M 


RNor 


15323- 


-4920 


M 


T Nor 


15402- 


-5449 


M 


UHer 


16235+1900 


M 


RR Sco 


16534- 


-3030 


M 


RT Sco 


17001- 


-3651 


M 


RAql 


19039+0809 


M 


TPav 


19451- 


-7153 


M 


SPav 


19510- 


-5919 


M 


RR Sgr 


19528- 


-2919 


SR 


RUSgr 


19552- 


-4159 


M 


RTSgr 


20144- 


-3916 


M 


UMic 


20259- 


-4035 


M 


R Aqr 


23412- 


-1533 


M 



363 


6.7 


315 


0.88 


332 


3.6 


128+18 


0.58 


235 


7.5 


592 


0.61 


408 


6.0 


308+102 


0.18 


278 


7.6 


593 


0.01 


391 


7.9 


493 


0.72 


171 


7.1 


532 


0.20 


427 


6.8 


250+53 


0.22 


368 


6.3 


238 


0.01 


341 


7.0 


351 


0.15 


141 


2.6 


61+5 


0.94 


309 


4.6 


128+14 


0.12 


150 


5.7 


255 


0.77 


531 


10.0 


380 


0.37: 


315 


8.0 


613 


0.33 


220 


7.0 


418 


0.94 


308 


7.9 


559 


0.92 


389 


4.5 


108 


0.63 


361 


6.0 


115+14 


0.44 


546 


6.3 


172 


0.01P 


218 


7.5 


570 


0.12 


508 


7.2 


195+71 


0.82P 


241 


7.4 


454 


0.33 


406 


7.5 


399 


0.35 


281 


5.9 


259 


0.55 


449 


8.2 


433 


0.11 


284 


6.1 


211+53 


0.43 


244 


8.0 


708 


0.51 


381 


7.3 


378 


0.82 


336 


6.8 


348 


0.69 


241 


7.2 


492 


0.82 


306 


7.0 


133+59 


0.86 


334 


8.8 


882 


0.33 


387 


6.5 


197+122 


0.64 



0.26 


0.25 


13,100 


21 


0.66 


0.25 


2,200 


7 


0.30 


0.48 


89000 


118 


0.27 


0.28 


14000 


22 


0.28 


0.25 


46,500 


64 


0.30 


0.30 


38,600 


54 


0.28 


0.25 


37,400 


52 


0.22 


0.19 


6,300 


13 


0.23 


0.18 


5,400 


11 


0.36 


0.30 


19,600 


30 


0.28 


0.26 


500 




0.26 


0.24 


2,100 


7 


0.24 


0.30 


10,300 


18 


0.21 


0.21 


16000 


25 


0.19 


0.22 


43,700 


60 


0.22 


0.24 


22,200 


33 


0.96 


0.33 


54,500 


74 


0.22 


0.21 


1,300 


6 


0.22 


0.23 


1600 


7 


0.63 


0.27 


4,200 


10 


0.18 


0.18 


30,900 


44 


0.22 


0.20 


4000 


10 


0.27 


0.24 


26,200 


38 


0.16 


0.17 


14,300 


23 


0.25 


0.24 


8,500 


15 


0.27 


0.24 


23,800 


35 


0.19 


0.17 


4000 


10 


0.25 


0.25 


66,300 


89 


0.24 


0.21 


15,900 


25 


0.25 


0.23 


14,700 


23 


0.21 


0.20 


25,600 


37 


0.23 


0.22 


2,100 


7 


0.20 


0.19 


78,200 


104 



detected (see text) 



*Distances with formal errors have been determined using Hipparcos trigonometric parallax measurements. Other distances have been 
determined from period-luminosity relations. 

"Brightness temperatures are calculated for an adopte d radio size of 5 X 10 13 cm (Section 4.2). 
f Shock velocities estimated using the model o f lFadevev & GilleJ Jl998l) . 
p primary maximum 



source visibilities were Fourier tra nsformed using 'natural we ight- 
ing'. The images were CLEANed ( lHogbomlll974|-|Clarkll980h and 
restored with a synthesized beam of minimum FWHM 1 arcsec at 
3 cm and 2 arcsec at 6 cm. The rms noise levels in the images were 
measured empirically from the images. 



Radio continuum emission was detected from one source only, 
the symbiotic binary R Aqr. This source is discussed below. For the 
remaining 33 sources, columns 8 and 9 of Table 1 list the 3-a upper 
limits to the 4.80 and 8.64 GHz flux densities which are between 
0.16 and 0.96 mjy at 4.80 GHz and between 0.17 and 0.48 mjy at 
8.64 GHz. 



4 DISCUSSION 



4.1 R Aquarii 

R Aquarii is a nearby symbiotic binary system with a cool Mira 
star and a hot dwarf or sub-dwarf companion star. The latter is evi- 
dent from strongly ionised optical emission lines but is not directly 
visible. The separation of the two stars has been estimated to be ~ 
15 AU and the temperatures of the Mi ra and companion stars to be 
~ 2 8 00 and 40000 K respectively ( iBurgarella. Vogel & Parescd 
1 19921 ; iHollis. Pedeltv & Lvonl [199% . Mass transfer between the 
two stars probably occurs as the hot companion captures mate- 
rial from the cool circumstellar wind of the Mira variable. Maser 
emission from SiO, H2O and OH molecules clearly indicates the 
presence of cool dusty ma t erial associated with mass-lo ss from 
the Mira dZuckermanl 1 19791 ; ISeaouist. Ivison & Halllll995l) . while 
a complex ionised envelope, detected in both radio continuum and 
optical emission, extends to at least 15 000 AU from the binary pair. 
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Figure 1. Radio continuum image of R Aqr at 6 cm, from data taken with 
the Australia Telescope Compact Array in November 1995. The contour 
levels are at -0.24, 0.24, 0.48, 0.96, 1.92, 3.84, 5.76 and 7.68 mJy/beam. 



Figure 2. Radio continuum image of R Aqr at 6 cm, from data taken with 
the Australia Telescope Compact Array in November 1995. The contour 
levels are at -0.24, 0.24, 0.48, 0.96, 1.92 and 3.84 mJy/beam. 



R Aqr also contai ns an astrophysical jet, first dete cted from 
optical spectroscopy by ( Wallerstein & Greenstein 1980) and in ra- 
dio continuum by dSopka et alJI 19821) . HST images have revealed 
a narrow, collimated and curved jet structure which extends from 
with in the orbit of the bina ry pair to a distance of several thousand 
AU jParesce & Hacklll994|) . 

Radio continuum emission was first detected from R Aqr i n 
1971 with a flux density of ~ 200 mjy at 85 GHz dJohnsod[l98(il) . 
iGregorv & Seaquistl dl974t) monitored the source during 1973 and 
found the radio emission to be strongly variable, on a timescale of 
months. The available radio flux measurements of R Aqr (includ- 
ing our ATCA result) indicate that the radio continu um emission 
has be en in a quiescent state since the mid-1970's. Iwillson et al.l 
(1981) suggested that the initial variability may have been due to 
the periastron of the companion star to R Aqr, revolving around the 
Mira in a highly eccentric orbit with P or b — 44 years. This result is 
suppo rted by VLA 43-GHz observations of lHollis. Pedeltv & Lvorj 
dl997l) . They spatially resolved the Mira variable (detected from 
its SiO maser emission) and the hot companion (detected in ra- 
dio continuum) and determined the spatial offset between the two 
stars to be 55±2 milliarcsec. Based on an orbital period of 44 
years, Hollis et al. estimate the semi-major axis of the orbit to be 
~ 2.6 x 10 14 cm with a high eccentricity of e = 0.8. 

The radio continuum emission from R Aqr has 
previously been imaged using several diffe rent config- 
urations of the Very Larg e Ar ray dSopkaetal] 1 19821; 



Kafatos, Hollis & Michalitsiano I 1 19831; | Kafatos et al.l 



19891; iHollisj Il985l; IHollis. Dorband & Yusef-Zadehl 
Hollis ^Pedeltv & Lvonl 1 19971 and references therein) 



1992; 



On 



the largest scale, the 6- c m im ages of IHollis et alj dl987l) : 
IHollis, Oliversen & Wagner! dl989h show optically-thin, thermal 
emission over an extended 2-arcmin region which is associated 
with the optical emission nebula. On much smaller angular scales, 
radio images at 2 and 6 cm show discrete components corre- 
sponding to compact regions within the jet. The strongest radio 
continuum emission occurs from a region located close to the Mira 
position, which has a spectral index of a — 0.5 (where S v oc v a ), 
consistent with optically-thick thermal free-free emission from 
an outflowing wind. Other jet components however have negative 
radio spectral indices indica ting non-thermal emission. 

iDoughertv etal] d 19951) observed R Aqr on MERLIN at 5 GHz 
(December 1992) and 1.7 GHz (August 1993). In their 5-GHz im- 
ages of R Aqr (resolution 40 milliarcsec), the jet is resolved into 
six, apparently thermal components with a total flux density of 
16.1±0.6mJy. 

In 1995 November the total 6-cm flux density detected with 
the ATCA was 17.2 mJy (Figure [T), in good agreement with the 
MERLIN value. The total detected 3 cm flux density was 15.9 mjy 
(Figure [2}, corresponding to a spectral index, averaged over the 
detected region, of —0.13. The ATCA radio images show the radio 
continuum emission from R Aqr to be slightly extended in a north- 
south direction. A comparison of the ATCA images with the far 
more sensitive earlier VLA and MERLIN images shows that radio 
continuum was detected from the inner jet region. The ire-coverage 
and sensitivity of the ATCA survey data were however not sufficient 
to resolve the complex structures within the jet, or to detect the 
extended envelope around the binary system. 

We do not discuss R Aqr further in this paper. 
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4.2 Brightness Temperature Limits 

For optically-thick radio continuum emission at centimetre wave- 
lengths, the gas brightness temperature, Tb, is related to the radio 
flux density, S v , by the Ray leigh- Jeans approximation so that 



S v = 



kT b , 



(3) 



where r s is the radius of the radio-emitting region, v is the fre- 
quency of the radio emission, D is the source distance and c is the 
speed of light. This may be written as 



T b = 3.95 x l(f 



D 



100 pc 



5 x 10 13 cm 



(ghs) 



(4) 



(mJy) K ' 



Direct optical interferometric measurements of AGB stars 
have determined photospheric radii of typically (2-3) xlO 13 cm 
dTuthill et al.ll 19941 : lHaniff, Scholz & Tuthilllll995h . while thermal 
continuum emission is expected to occur from outside the optical 
photo sphere but within the grain condensation radi i of a few stellar 
radii i Greenhill et alj[l995l) . From their radio data lReid & Mentenl 
i ll 9971) ' have shown that at cm-wavelengths the thermal radio emis- 
sion from an extended photosphere reaches unity optical depth at 
a radius of approximately two stellar radii. For comparison with 
their results, column 10 of Table 1 gives upper limits to the bright- 
ness temperatures determined from equation (4) using the stellar 
distances given in column 6 of Table 1 and a 'characteristic' radio 
size of r s 



4.3 Shock Velocity Limits 

For optically-thick thermal radio continuum emission, we can use 
the observed brightness temperature, Tb, as a direct measure of the 
post-shock electron temperature T e . To obtain 'reasonable' esti- 
mates of the shock velocities we adopt a characteristic size for the 
radio-emitting region of 5 x 10 13 cm. 

Column 1 1 of Table 1 gives the upper limits to the shock ve- 
locities determined using the model of ( Fade vev & Gilletl 19981) and 
up per limits to the radio bri ghtness temperatures given in column 
10. IFadevev & Gillel dl998l) obtained a self-consistent solution of 
the equations of fluid dynamics, radiation transfer and atomic ki- 
netics, involving parti al ionisation and excitati on of bound atomic 
states. In Table 2 o f IFadevev & Gilletl dl998h the dependence of 
the post-shock electron temperature T e on the shock velocity v s is 
nearly linear and can be approximated as 



778v s - 3337. 



(5) 



Shock velocity limits have been calculated using formula Q. 

The primary result of our ATCA survey was the lack of any 
radio continuum emission detected from the 33 sources observed 
(excluding the symbiotic star R Aqr). The brightness temperature 
limits given in Table 1 cover a wide range of values, due primarily 
to their strong dependence on the stellar distances. 

For the six nearest sources in our sample (L2 Pup, R Hya, 
W Hya, o Cet, R Car and RT Sgr), with distances within 140 pc, 
the upper limits to the brightness temperatures are comparable to or 
lower than the stellar photospheric temperatures. For these sources, 
we can exclude strong shock activity near two stellar radii with 
shock velocities above ~7 km s _1 . For the nearest source in our 
sample, the semi-regular variable L2 Pup, the lack of any radio 



continuum detection is surprising. For this source the upper limit 
to the 3-cm radio flux density is a factor of three lower than pre- 
dicted by the Reid & Menten's model and a factor of nearly two 
lower than expected for blackbody emission from an optical pho- 
tosphere of radius 3 x 10 13 cm and temperature 2 500 K. A likely 
explanation is that the stellar size is smaller than the adopted value 
used to estimate brightness temperatures. Our flux density limits 
are consistent with a stellar radius of 1.5 x 10 13 cm or smaller. For 
the other five sources with distances below 140 pc, the brightness 
te mperature limits are c onsistent with the radio-photosphere model 
of lReid& Mentenl dl997b . 

For a further five sources, R Lep, U Ori, R Cen, R Nor and 
R Aql, which are within 250 pc, the brightness temperature limits 
are below 6,300 K, consistent with stellar shock velocities below 
10-13 km s _1 . We note that these are conservative upper limits and 
that actual shock velocities may be considerably lower. 

For all 15 sources in Table 1 with distances within 320 pc, we 
can rule out the presence of strong stellar shocks with velocities 
above 22 km s _1 . For the 17 sources in the sample, with distances 
greater than 320 pc, our data are insufficiently sensitive to provide 
useful limits to the gas temperatures or shock velocities yield higher 
upper limits for us, ranging between 23 and 118 km s _1 . 

4.4 Comparison with Previous Radio Detections 

Many searches have been made for radio continuum emission 
from long peri od variable stars with largely negative results (e.g. 
IWendkenll995l and references therein). In Table 2 we summarise 
previously published detections of radio continuum emission at cm- 
wavelengths from Mira and semi-regular variable stars. 
The columns of Table 2 are: 

(1) source name (variability type, spectral type) 

(2) D, the stellar distance. For o Cet, W Hya, V Hya and R Aql, 
distances are as for Table 1. For CW Leo, the distance is taken from 
Crosas & Menten (1997). For RZ Ari the distance is determined 
from the Hipparcos parallax. Other distances are taken from the 
references given in column 8. 

(3) the year of observation 

(4) v, the frequency (GHz) 

(5) Sv, the detected flux density (mJy) 

(6) radio telescope used 

(7) reference for radio detection 

(8) Tb, the radio brightness temperature calculated assuming op- 
tically thick radio emission for a 'characteristic' radius of 5 x 
10 13 cm (Section 4.2) 

(9) v s , the shock v elocity calculated from the model of 
IFadevev & Gilletl dl998l) . using the brightness temperature given in 
column 8. 

Apart from R Aqr (Section 4.1), radio continuum emission has 
previously been detected from four of the sources in our sample, 
V Hya, R Aql, o Cet and W Hya. 

For the carbon star V Hya, our three-sigma upper limit of 
0.21 mJy is consistent with the 1989 VL A 8.4-GHz detection of 
0.22 mJy bv lLuttermoser & Brownl dl992l) . 

For R Aql, strong flare activity and non-thermal emis- 
sion has been reporte d from seve ral single-dish observations: 
IWoodsworth & Hughes! dl973lll977h observed a rapid rise in the 
10.6 GHz emission from 23 to 240 mJy over a period of only 
30 minutes. Less intense, but sti ll unusually strong emissi on from 
R Aql was later observ e d by iBowers &Kundj Jl979l) . while 
Estalella, Paredes & Rius ( 1983) detected non-thermal emission 
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Table 2. Mira and semi-regular variable stars with previous detections of radio continuum emissionat cm-wavelengths 



Star 


D* 


Year 


V 




Tel 


Ref 


n 


v s 


(Var, Sp) 


(pc) 




(GHz) 


(mJy) 
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(kms- 1 )! 
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-K 7 


744-0 9^ 
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8 4 
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1 900 
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Q 
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IS 4Q 
u. j J— u.*+7 
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2 
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7 
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Q44-0 1 1 
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Q 
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9 004 /OS 
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4 
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22 5 


i . i j 
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4 
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4 
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U. i o— u.z / 


VLA 
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7 
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1 47 ±0 20 
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6 
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j 


15 500 


24 






1 Q87 
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15 
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It. y 


9 1 44-0 1 8 2 
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Q 
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U. 1 ^tlCU.U J 


VLA 
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10 
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1 Q8Q 

1707 


8 4 


14 
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7 


1 000 

1 ,uuu 








1987 
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494-0 1 9 
U.M-Z^U. 1 z 


VLA 


7 


Q00 

7UU 


6 
O 


ci 1 Her (SR, M) 


117 ±38 


1983/84 


5 


0.83±0.06 2 


VLA 


5 


18,000 


28 






1984 


15 


1.71±0.10 


VLA 


5 


4,100 


10 


R Lyr (SR, M) 


107 ±6 


1987 


15 


0.45±0.12 


VLA 


7 


900 


6 


R Aql (M, M) 


211 ±53 


1982 


2.3 


381±16 


DSN 64 m 


11 


see 


text 






1982 


8.4 


40±6 


DSN 64 m 


11 










1970 


10.5 


23^240 


ARO 46 m 


12 
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14.9 


5.3±2.0 


Bonn 100 m 


13 










1989/90 


14.9 


0.24±0.18 


VLA 


3 


1,900 


7 






1989/90 


22.4 


0.80±0.4 


VLA 


3 


2,800 


21 


X Cyg (M, S) 


135 ±11 


1989/90 


8.4 


0.27±0.04 


VLA 


3 


2,800 


21 






1989/90 


14.9 


0.79±0.20 


VLA 


3 


2,600 


21 






1989/90 


22.4 


1.62±0.30 


VLA 


3 


2,300 


20 


R Cas (M,M) 


160 ± 13 


1989/90 


8.4 


0.11 ±0.05 


VLA 


3 


1,600 


7 






1989/90 


22.4 


0.79±0.30 


VLA 


3 
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7 



Notes to Table 2: 

Other names: a IRC+10216, AFGL 1381; 6 CIT 6, IRC+30219, AFGL 1403; 

1 Range of flux densities obtained from VLA monitoring data 

2 Average of two or more values 

3 Mira A only 

Re ferences: 

(lll Spergel, Giuliani & Knapd jl983l). (2)lDrake, Linsky & Judge! j 19931). (3) lReTd & Mentenl 4l997l>. (4)lMatthews & Karovskd i20O6h. 
(5) 



Drake & LinskvHl986l) (6)llCnapp et all il995l). (7) 



Drake et al 



199 lh . (8) Sahai. Claussen & Massonlh989t), (9)lMenten et all J200i 



(lO^ lLuttermoser & Browd <1992l) , (1 D lEstalella, Paredes & RiusHl983l) , (121 1 Woods worth & HughesHl973|) , (13^ lBowers & KunduHl97' 



L97g) . 
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with a spectral index, between 2.3 and 8.4 GHz, of —1.73. A pos- 
sible explanation for such a steep spectral index is non-thermal 
synchr oton emission fro m electrons accelerated by a strong stellar 
shock (Rudnitskij 1993). The more recent interferometric data for 
R Aql however shows no evidence for excess continuum emission. 
Our 8 GHz upper limit o f 0.2 mjy is consistent with the 1990 VLA 
upper limit of 0.09 mjy jReid & Menten|[l997h . The association of 
the earlier single-dish detections (which did not have accurate ra- 
dio positions) with the stellar source is unclear. We do not however 
rule out the possibility that unusual and irregular flare activity may 
occur from R Aql. 

For W Hya, our upper limit of 0.25 mjy i s slightly lower 
than th e mean flux density of 0.36 mjy obtained bv lReid & Mentenl 
( 1997). Their observations partially resolved the stellar disk, with 
an average diameter of 0.08 ±0.02 arcsec and a brightness temper- 
ature of 1,500 ± 570 K . 

For o Cet, the three-sigma detection of 

Sperg el. Giuliani & Knap^ dl983h_ is a marginal res ul t. As 
this is inconsistent with both Drake, Linskv & Judge! h993h 
and IReid & Menten] Jl997l), the published detection is probably 
spurious. IReid & Mentenl jl997l) found marginal evidence for 
periodic variability in the radio continuum emission with a mean 
flux density at 8.4 GHz of 0.42 mjy and peak-to-mean ratio of 
~ 1.1. From a periodic fit to the radio fluxes they showed that 
such variability could be explained by a modulation of the thermal 
radio emission caused by a weak pulsation-driven stellar shock 
propagating through the radio-photosphere, with a shock velocity 
of approximately 7 kms -1 . Our 3-sigma upper limit at 8.6 GHz 
of ~ 0.25 mjy is lower than the 8.4 GHz flux densities of 0.35 to 
0.49 mjy obtained by Reid & Menten in 1990-1991. The cause 
of the discrepancy between the VLA flux densities and the ATCA 
flux density limit, obtained five years later, is unclear. However, it 
is well-known that pulsations amplitudes and shock strengths vary 
from cycle-to-cycle and some variation in the radio continuum 
emission appears likely. Source variability may als o occur on a 
longer timescales. iTuthill. Haniff & Baldwin! i 19951) showed that 
the optical size of o Cet, measured at wavelengths between 700 
and 900 nm, increased systematically by a factor of 1.8 between 
1991 September and 1993 December. The increase in optical size 
may imply an episode of increased mass-loss and a build-up of 
the circumstellar dust shell. Such an enhancement of circumstellar 
dust density ma y damp the shoc k waves in the inner part of the 
stellar envelope (Rudnitskii 1 997n . 

The star o Cet is a symbiotic binary comprising a variable red 
giant (Mira A) and an companion (possibly a white dwarf, Mira B) 
accreting m aterial from the primary. It w as resolved in VLA obser- 
vations of Matth ews & Karovskal d2006h . The radio fluxes at 8.5, 
14.9, 22.5 and 43.3 GHz were measured separately for both com- 
ponents at several epochs in 2004—2005. In Table 2 we list only the 
fluxes measured for the red-giant component Mira A. 

For 1 1 of the sources listed in Table 2 (o Cet, p Per, /i Gem, R 
Leo, RW LMi, W Hya, g Her, R Lyr, R Aql, \ Cyg and R Cas) the 
detected radio continu um flux densities agre e well with the radio- 
photosphere model o fReid& Menten (1997). For these sources we 
estimate brightness temperatures of 1500— >4000 K indicating little 
if any excess emission. From the shock model given in paper I, 
the radio brightness temperatures indicate that only weak or mod- 
erate shock strengths are possible with shock velocities below ~ 
lOkms- 1 . 

As discussed by lReid& Mentenl dl997t) . at temperatures be- 
tween ~ 1 000 and 4 000 K, free electrons occur primarily from 
elements such as potassium, sodium and calcium which, although 



far less abundant, have much lower ionisation potentials than hy- 
drogen. Their calculations show that the dominant source of the 
free-free opacity occurs from interactions between free electrons 
from ionised metals with neutral hydrogen atoms and molecules. 

4.5 AGB Stars with Excess Radio Continuum Emission 

Excluding R Aql and o Cet (discussed above), only three sources 
in Table 2, CW Leo, V Hya and a 1 Her show a strong excess of ra- 
dio continuum emission above the levels expected for radio photo- 
spheres. We note that two of the three sources, CW Leo and V Hya 
are N-type AGB carbon stars. The carbon stars represent around 
eight per cent of AGB stars and have C/O ratios enhanced surface 
carbon abundances due to a dredge-up of processed material. These 
stars have high mass loss rates and are considered likely to be near 
the end of the AGB stage of evolution. 

CW Leo (=IRC+10216) is an extreme carbon star (C9) and 
a strong source of molecular emission, with over 50 molecular 
species detected from i ts circumstellar envelope . From an analy- 
sis of the CO properties. ICrosas & Mentenl (l997) estimate the dis- 
tance to be 150 pc. For this distance, the cm radio continuum flux 
densities indicate radio brightness temperatures between ~5,000 
and 15,000 K. For our model, the corresponding shock velocities 
are between 11 and 24 km s , For CW Leo, the radio continuum 
emission almost certainly occurs from a chromospheric region. IUE 
ultraviolet spectra with ionised emission features from Mgll, CII 
and Fell have shown that both ox ygen-rich and carbon-rich AGB 
stars have stellar chrom ospheres jjohnson & Luttermoseil Il987l) . 
Luttermoser et al. Ul989h provided a chromospheric model for re- 
type carbon stars. In their model the inner chromospheres have a 
steep temperature gradient and an outwards expansion velocity of 
around 50 km s _1 , while the outer chromospheres are nearly static 
relative to the stellar photospheres. The heating mechanisms in the 
chromospheres of late-type stars are not well understood. However, 
if strong shocks are present then some chromospheric heating is 
likely to occur due to th e dissipation of sh ock energy and heat- 
ing of the post-shock gas. lMenten et alj d2006h monitored CW Leo 
in 1991-1993 in the radio continuum with VLA at 8.4, 14.9 and 
22.5 GHz (see Table 2). They found possible variability correlated 
with the infrared phase and a spectral index very close to 2. The 
variability, observed flux densities and upper limit on the source 
size are consistent with the emission arising from the stellar pho- 
tosphere or a slightly la rger radio photosphere. Using the model of 
iFadevev & Gilletl dl998l) . we have estimated shock velocities that 
would heat the radio photosphere to temperatures listed in Table 2; 
we have obtained v B = 14-15 km s _1 . 

V Hya is also an extreme (C9) carbon star. CO imaging of 
V Hya has shown that most of the stellar ma ss-loss occurs withi n 
a high- velocity wide-angle bipolar outflow dKahane et alj|l99fib . 
Such a geometry indicates that the star is likely to be close to leav- 
ing the AGB and may be in a short-lived common envelope binary 
system. iKahane et alj d 19961) find that along the bipolar axis there 
is a strong velocity gradient with velocities decreasing outwards 
from over 50 kms -1 at the inner boundary to around 8 kms -1 
at the outer boundary. Within the bipolar cones, multiple shocks 
are likely to occur as faster moving gas collides with more slowly 
moving material. From this model and the above discussion we 
sugge st that the radio continuum emission from V Hya, detected 
by iLuttermoser & Brownl dl992l) . has two possible sources; chro- 
mospheric emission from the outer atmosphere and emission from 
high-velocity shocked regions located near the base of the bipolar 
outflows. 



8 Georgij M. Rudnitskij and Jessica M. Chapman 



Finally, a 1 Her is an M5II star which belongs to a vi- 
sual binary system. For this source the 5 — >15 GHz spectral in- 
dex of 0.8 is close to the canonical value of 0.6 for thermal 
emis sion from an ionised w i nd expanding at a cons tant velocity 
(e.g. lWright & Barlowlll975b . lDrake & Linskvl ( fl986b interpret the 
emission from this source as thermal emission occuring from a 
chromospheric region of temperature ~ 10,000 K and radius ~ two 
stellar radii. Here we note that the spectral index of 0.8 is consistent 
with expansion in the stellar chromosphere. 



4.6 Comparison of Radio and Optical Studies 

Considerable controversy exists on the interpretation of optical 
spectra of M-type Mira and semi-regular variables. These char- 
acteristically show hydrogen Balmer and metallic line emission 
lines, which are strongest near maximum light and decline in inten- 
sity with increasing phase towards minimum light. The linewidths 
of the emission lines are also strongly phase dependent with full 
widths of typically 80 km s -1 near light maxim u m and 40 km s -1 
near light minimum (e.g. iGillet et all Il985 altl iRichter & Wood 
1200 lb . 

From an analysis o f the Balmer linewid t hs, Gillet et al.l 

jl985albT) ; lGliieil<1988allbh : lGillet. Lafon & David < fl989h advocate 
a model in which the double-peaked structure of the Ha emission 
line profiles observed in some M-type Miras occurs from the front 
and back-side of a single spherical shock wave, with shock veloci- 
ties es timated, in most cas es, to be between 40 and 70 kms -1 . For 
S Car. lGillet eta l. (1985b) estimated a maximum shock velocity of 
90 kms" 1 . 

Other authors determine smaller shock velocities by argu- 
ing that multiple shocks exist simultaneously in Mira atmospheres 
(Wood 1979l lBoweiJl988l : lDerviz & Somsiko\ll99Cl : IWoodsworthl 
1995). In thi s case, the Ha emis sion line profiles may have com- 
plex shapes. IWoodsworthl |l995) found three components in the 
Ha line profiles of several S-type Miras. These were modelled as 
emission from different ionised layers, associated with a 'main' 
outward shock and an 'inner' shock caused by infalling post-main- 
shock gas. The shock velocities, inferred by this model, are below 
40 kms -1 . 

From the radial velocities of optical h ydro gen and metallic 
emissi on lines, iFox. W ood & Dopita ( 1984) and IRichter & Woo3 
j200 lh showed that at maximum light the post-shock emission re- 
gion has an outwards velocity of 10-20 kms -1 . Similar veloci- 
ties h ave also been obtained from low-opacity ultr a violet emission 
lines dWood & Karovskal l2OO0b . IRichter & Woo3 j200lh demon- 
strate that the radial velocities decrease with increasing phase, and 
are almost zero at minimum light. This can be interpreted as a de- 
crease in shock strength as a shock wave moves outwards through 
the stellar atmosphere, while the decrease in linewidth is likely to 
reflect a decrease in temperature and hence Doppler broadening be- 
hind the shock front. 

In the present study we have argued that the lack of any de- 
tected radio continuum emission from the 1 1 nearest sources in our 
sample is consistent with an upper limit to stellar shock velocities 
of around 6-7 kms -1 near two stellar radii. A similar conclusion, 
based on the lack of s trong radio variability, was also reached by 
lReid& Mentor] ( Il997h . The lack of radio continuum detections in 
the present study clearly favours low shock velocities and/or damp- 
ing of the shock energy within the inner few stellar radii. 



5 CONCLUSIONS 

We have made a sensitive search with the ATCA for radio contin- 
uum emission at 3 and 6 cm from 34 Mira and semi-regular variable 
stars. Our search provided no detections of radio continuum emis- 
sion, apart from for R Aqr which is a well Known symbiotic binary 
system. In 1995 November the radio emission from R Aqr was in a 
quiescent state. 

From the upper limits to the flux densities determined for the 
other 33 sources, we have estimated upper limits to the gas bright- 
ness temperatures near two stellar radii at a characteristic size of 
5 x 10 13 cm. Shock velocities have be en estimated using the model 
presented bv lFadevev & Gilled dl998h . 

For the six nearest sources in the ATCA survey, with distances 
within 140 pc, we obtain brightness temperatures below 2,300 K 
indicating that there is no excess of radio continuu m above the lev- 
els exp ected from the radio-photosphere model of iReid & Mentenl 
( 1997). From previously published detections of centimetre contin- 
uum emission we find that for 10 out of 13 sources the detected flux 
densities indi cate radio brightness tem peratures below 4,000 K. In 
the model of iFadevev & Gilled Jl998b we have used, a brightness 
temperature limit of 4,000 K corresponds to a shock velocity limit 
of ~10 kms -1 . 

Only three long-period variable stars of those observed by 
other authors, CW Leo, V Hya and a 1 Her, show a strong excess 
of radio continuum emission. In each case, the strong continuum 
emission is likely to be associated with a stellar chromosphere. 
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